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Abstract 

Preliminary results toward the analysis of the Hamiltonian structure 
of multifield theories describing complex materials are mustered: we in- 
volve the invariance under the action of a general Lie group of the balance 
■ of substructural interactions. Poisson brackets are also introduced in the 

material representation to account for general material substructures. A 
Hamilton- Jacobi equation suitable for multifield models is presented. Fi- 
nally, a spatial version of all these topics is discussed without making use 
of the notion of paragon setting. 



^ ; 1 Lagrangian and Hamiltonian descriptions of 

elastic bodies with substructure 

In standard continuum mechanics, each material element of a body is "col- 
lapsed" into the place occupied by its centre of mass; let X be that place in 
the reference placement; the set of all X is taken to be a fit region Bq of the 
Euclidean space. 

Sometimes such simplicistic model of physical reality is insufficient; then, to 
render the picture adequate, the material element must be portrayed as a system 
and at least some coarse grained descriptor v (an order parameter) enters the 
picture. 

Here, as in [1] (see for other details and additional results [2-5]), we take 
1/ as an element of some differentiable manifold M , and presume that physical 
circumstances impose a single choice of metric and of connection for A4 . 
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We also assume that the region occupied by the body in the current place- 
ment be obtained through a sufficiently smooth mapping x -.Bq ^ £ ; so that the 
current place of a material element at X in Bq is given by x (X); and B = k (Bq) 
is also fit; we denote as usual with F the placement gradient. We presume also 
that another sufficiently smooth mapping u iBq — > A4 shows the present value 
of the order parameter at X. A motion is a pair of time-parametrized fami- 
lies Xt (X) = x(X, and Vt (X) = i>'(X,t), twice differentiable with respect 
to time. Rates in the material representation are indicated with x (X, t) and 
i>(X,i). 

We restrict here our attention to bodies for which a Lagrangian density C 
exists, so that the total Lagrangian L of the body is given by 

Lb„ - / £(X,x,x,F,iv,«>,Vi^). (1) 

"'Bo 

the covariant gradient Vi/ being based on the mandatory connection. We pre- 
sume that C be of the form 

1 2 

C (X, X, X, F, v, i>,\/v>) = ||x|| + («^, i') ~ Po^ (X, F, u.Vu) - p^w (x, v) , 

(2) 

where pg is the referential mass density (conserved during the motion), x the 
kinetic co-energy (see [1], p. 19) associated with the substructure, e the elastic 
energy density and w the density of the potential of external actions, all per unit 
mass. Below we use the notation b — —dy^w for the density standard external 
actions and (3 — —diyW for the substructural ones. The kinetic energy density 
PqK (u, i>) pertaining to the substructure is the partial Legendre transform of x 
with respect to v>. 

If C is sufficiently smooth, we may apply standard procedures to derive 
Lagrange equations for the functional L^o'- 

d^^d^C- DivdvC (3) 



doC = d^C - DivdvuC. (4) 

where Div is the divergence calculated with respect to X, i.e. Div = trV. 
Put 

— x-dxjC + i/doC — C. (5) 

Clearly, is the density of the total energy. In fact, since di,C — Pi^doX: ttiG 
term i> ■ doX ~ X ^ coincides with the substructural kinetic energy density 
K, (f , i>) (hence the presence of x rather than k, in the expression of C), then 

1 2 

^l=2^ol|x|| +poK(i/,i>)+Poe(X,F,j^,Vi^)H-pow(x,jy), (6) 
as asserted. 
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The balance of energy can be expressed in terms of Sj as follows 

h - Div (xP + i^S) = 0, (7) 

where P and S are respectively the Piola-Kirchhoff stress and the referential 
microstress 

P = -dFC , S = -dv.C. (8) 
That lO is true follows from direct computation. 

Equations Q and l@J lead us to an appropriate version of Nother theorem 
(see [4] p. 29); here we follow the program of [6] (p. 284). We consider some 
virtual motion of our system, by assigning two one-parameter families f]. of 
sufficiently smooth point valued diffeomorphisms, i = 1,2, acting respectively 
on Bo and £, and a Lie group G of transformations of A4. We indicate with a 
prime the derivative with respect to the relevant s. 

1. At each si, f}^ acts on Bq so that X i — > f^^ (X) £ £, and is isocoric (no 
virtual change of density), i.e. Divil[ = 0; fp is the identity. We put 
fo^' (X) = w (X), then w'= 0. 

2. At each S2, f^^ is a diffeomorphism that transforms £ into itself. We 
assume that fg is the identity and put fp ' (x) ~ v (x) . 

3. A Lie group G, containing SO (3), acts on M and the infinitesimal gener- 
ator of this action at v is indicated with {v) (see [7], p. 256); Ug is the 
value of V after the action of 5 G G. If we consider a one-parameter tra- 
jectory S3 G G such that go is the identity, we have also S3 ^ ^qs^ 
and {u) = ^I'g^g |s3=o- When G coincides with the special orthogo- 
nal group SO (3), we identify {u) with ylq, being q the characteristic 
vector of a rotational rigid velocity and A a linear operator mapping vec- 
tors into elements of the tangent space of M , namely, if I'q is the value 
of the order parameter measured by an observer after a rotation q, then 

Henceforth, to simplify notations, we use f^, and Vg to indicate f,^^^ (X), 
(x), i/g^, (X), and write |o for |si=o,s2=o,s3=o • Moreover, grad indicates the 
gradient with respect to x. 

We say that C is invariant with respect to P's and G when 

£(X,x,x,F,«^,i>,V«^) = 



-£(f^f^(5rao^f')x,(5radf')F(Vfl) \ i^g, J>g, (Vi.,) (Vf i) . (9) 
Let us define 

Q = diC- (v - Fw) + dcC- i'') - (Vi^) w) (10) 
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^ = /:w+ [d^cf (v - Fw) + [d^^cf i^) - (vi^) w) , (11) 

where v, w and [u) are as mentioned in items 1, 2, 3. 

Theorem 1 (Nother-like theorem for complex materials). If the Lagrangian 
density C is invariant under f}_^ , f^^ and G, then 

Q + Divd - 0. (12) 
Proof. To prove the theorem, as a first step we note that Q implies 

(f\ f2, (5radf2) X, (5radf2) F (Vf^)"' , i^^, v,, (Viv^) (Vfi)-') |o = 0, 

(13) 

/:(f\f2,(5radf2)x,(5radf2)F(Vfi)"\i.,,i/,,(Vi.,)(Vfi)"') |o = 0, 



ds-2. 



(14) 



(f\ f2, (5radf2) X, (5radf2) F (Vf^) ' , iv^, (Viv,) (Vf^) |o = 0, 

(15) 

which lead to 

dy^H ■ w - dwC- (FVw) -dx^^C- (Vi^Vw) =0, (16) 



5x/: • V + 9x/:- ((57-adv) x) +dFC- {{gradv) F) =0, (17) 

d^C ■ {^) + doC ■ (i^) + d^^C ■ (i^) = 0, (18) 

as a consequence of the properties listed under items 1, 2 and 3 above. Then, 
we calculate the time rate of the scalar Q, the divergence of the vector ^ and, 
by using the equations (O and identifying S3 with t, we recognize that 

q + Div^^^C\q+^C\o+^C\o, (19) 
dsi ds2 dss 

which proves the theorem. □ 

Remark 1. As a first special case, we require that alone acts on C leaving 
V arbitrary. By using p7|) we obtain from p2l) 

\ 

g-^d^C - d^C + DivdwC j • V = 0, (20) 
i.e. 

PqX ^Poh+DivP, (21) 
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which is the standard equation of balance of momentum. 

Remark 2. With G arbitrary, we consider its action alone on £; by using 
p8|) . with the identification S3 — t, we obtain from (|12|l that 

(^§-^doC - d^C + Divdv^C^ ■ {^) = (22) 

or 

Po - duX^ + z-pof3- DivS = 0. (23) 

z — —p^duC is called self-force in the terminology of [1]. This result assures 
the covariance of the balance of substructural interactions. When G coincides 
with 5*0 (3), the co- vector into parentheses in H22(l . namely the term multiplying 
{iy>), must be an element of the null space of (see for details of this special 
case [1-4]). 

Remark 3. As a second special choice, let be such that v = q x (x — Xq) 
(with q a rigid rotational velocity xq a fixed point in space) and G — SO (3); 
qx is an element of its Lie algebra, thus (i^*) = .4q. If £ is independent of x 
and we assume that only f ^ and G (in the form just defined) act on £, we have 

skw {dpCF'^) = e (^A^d^C + {VA^Y d^^C^ , (24) 

where e is Ricci's alternating tensor and skw {■) extracts the skew-symmetric 
part of its argument. 

Remark 4. If we require that alone acts on C, with w arbitrary (but 
satisfying 1), by using (|16|l we obtain from (|12|l that 

{F^diL + Vi^^doC) - Div Qpo l|xf + X «>)) l) - d^C = 0. (25) 

where P = el — F-'^P — Vu^±S is the modified Eshelby tensor for continua with 
substructure (see [4] for a similar result in a non-conservative setting, where the 
elastic potential e is substituted by the free energy). I is the second-order unit 
tensor and the product * is defined by {Vu^ ±S) n • u = Sn- {Vu) u for any pair 
of vectors n and u. 

Remark 5. Let us assume as special choices that is such that w =q x 
(X — Xo) (with q a rigid rotational velocity Xq a fixed point in space) and that 
G — SO (3), being qx an element of its Lie algebra, thus {v) — Ai\. If the 
material is homogeneous, and we assume that f ^ and G alone (in the form just 

defined) act on £, we have skw {f'^DfC + {Vvf *d^^c}j = 0. 

Remark 6. The action of f ^can be interpreted as a special virtual mutation 
of a possibly existing smooth distribution of inhomogeneities throughout the 
body, in the sense of [8]. In other words, we may say that 1)25(1 is the balance of 
interactions arising when the body mutates its inhomogeneous structure. This 
interpretation has been also suggested in [9] in non-conservative setting. 
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1.1 Hamilton equations 

Define p and fi, respectively the canonical momentum and the canonical sub- 
structural momentum, by 

p =diC , n =dcC. (26) 

The Hamiltonian density 7i, 

n (X, X, p, F, V, ti,\/iy) = p-x + /X • z> - £ (X, X, x, F, v,Vv) . (27) 

has partial derivatives with respect to its entries; some of them are the oppo- 
site of the corresponding derivatives of C so that Q can be also written 
respectively as 

p = —dxH + DivdrH, 

X = dpH; (28) 
fi = —di^H + DivdxjvH, 

= d^H. (29) 

2 Canonical Poisson brackets in multifield the- 
ories 

We now consider a general boundary value problem where the following bound- 
ary conditions are associated with and 

x(X)=x ona^^^iSo, (30) 

dpHn = t on a'^'^o, (31) 

i/(X) = P ond'^^^Bo, (32) 

d^^Hn on d'^^^Bo] (33) 

X, t, D and t are prescribed on the relevant parts of the boundary and CI (OBq) = 
CI {d^^^Bo U 5(2)Bo), with d^^'^Bo n d'-^^Bo = 0, where CI indicates closure and 
n is the outward unit normal to dBo at all points in which it is well defined. 
We assume that there exist two surface densities U (x) and U (i/) such that 

t = Pod^U , t=Pf,d^U, (34) 

where U and U plays here the role of surface potentials. 
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Then the Hamiltonian H of the whole body is given by 
i/(X,x,p,«^,/x) - / H(X,x,p,«^,m)- / (C7 (x) - [/ (iv)) . (35) 

Notice that we write H (X, x, p, v>, fi) instead of H (X, x, p, F, v>, /^jV^") because 
below we consider directly variational derivatives. 
Theorem 2. The canonical Hamilton equation 

F^{F,H} (36) 

is equivalent to the Hamiltonian system of balance equations (I28|l - (|29|) for a con- 
tinuum with substructure where F is any functional of the type f (X, x, p, f , /x), 
with f a sufficiently smooth scalar density, and the Poisson bracket {•, •} for a 
complex material is given by 



{F,H} 



Bo 



6f sn sn Sf 

Sx Sp (5x Sp 

(5f sn I sn Sf . 

g(2)Bo \Sx dp ' dx dp ' 

rsi sn sn si 

\Sv Sfi Sfi Sv 

f (Sf sn, sn Sf , \ 

where the variational derivative ^ is obtained fixing p and allowing x to vary^ ; 
an analogous meaning is valid for the variational derivative with respect to the 
order parameter. 

The proof can be developed by direct calculation. Clearly, {■,■} is bilin- 
ear and skew-symmetric, and one can check easily that it satisfies the Jacobi's 
identity. We note that 

^ f -^■'«"«-™)- 

f (Sf dn Sf \ 
+ / ^ • — ^ ■ y"''^ ^ Divd^^n) + 

J Bo \St^ dfl Sfl j 

^ f fSf dn, (f, rr a -u ^\ ^ 

+ [s^ ■ dJI " sjl ■ (9.u~d^^nn) j 



^See relevant remarks in [10]. 
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and, in terms of functional partial derivatives, 

^ ffSf.Sf.5f. 5f . 

Js„ \STi Sp Su Sfi 

f Sf : f Sf : 

By identifying analogous terms in H38|l and H39|l . we obtain both the Hamiltonian 
system (PS)l - (|^ and the boundary conditions 

When we put F = H, H36() coincides with the equation of conservation of 
energy. We have, in fact, 

H = [H, H} = 0. (40) 

Geometrical properties of the Poisson brackets for direct models of rods, 
plates and complex fluids have been discussed in [10], [11]. 

3 A formal approach toward an Hamilton- Jacoby 
theory with gradient effects 

Let h he a, smooth diffeomorphism 

h: (X,x,p,F,iy,/x,Viv) ^ (X,x,,p,,F,,zy,,/x,,Vi/*). (41) 

The transformation h generates a new Hamiltonian density 

W4X,x„p,,F,,iy„/x„Vz/,), (42) 

with corresponding Lagrangian density 

= p* -x, +//^ • i>, -Ti:,. (43) 

If h were such that = 0, then an immediate integration of the system (|28|l . 
H29|l could be achieved. To this aim we choose h to be such that the integral of 
the difference C — C<, between two instants, say ti and t2, be equal to the time 
derivative of a generating function S of the type S — S {t, X, x, p^, i/, /x^), i.e. 

\-C,^S\t^t,-S\t^t,- (44) 



Then, from H44() we would have 

(p • X + /.t ■ i> - 7Y) - (p* • X* + /.t^ • - 7Y*) = 

= S ^ dtS + d^S ■ ± + dp,S ■ p^ + d^S ■ + d^^ S-fi,, (45) 

and hence 

p=d^S , tJ,=d^S, (46) 
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X* - xo = S* , V = d^^ S, 



(47) 



dtS + n = n^. (48) 

To obtain 147|l one makes use of the fact that S J^^ (p- (x — xq) + ■ v) — 
for variations vanishing at ti and t2 (in the sense that 5 (p ■ u) \ t=ti.t2 = ^-nd. 
S (p- (x — Xq)) \t=ti,t2 = 0) so that p ■ x = p- (x — xq) and fi ■ C> ~ fj, ■ v. 

A necessary and sufficient condition to assure that = is 

dtS + n (X, x, d^S, F, u, d^S,Vu) = 0, (49) 

which is a Hamiltonian-Jacobi hke equation. Since — 0, p^,, and fi^ are 
constant in time, the time derivative of S reduces to 

S = dtS + d^S ■ ±+duS ■ u = -n+p ■± + (j, u=C. (50) 

The relation H50() allows us to determine S to within a constant, namely 

S = Cdt + const. (51) 



4 The spatial form 

Circumstances in which the notion of reference placement is wanting, as in 
the case of fluids or granular flows, render the choice of a material or spatial 
representation not matter of form only (see, e.g., [12] for standard bodies). Here, 
having in mind the study of complex fluids, we provide a spatial variational 
derivation of the balance equations free of any concept of reference place or 
paragon setting and without even formal recourse to an inverse motion. So, 
in the present section x G B is just a point in space. The notation u (x, t) 
is used for the velocity field over B. The order parameter is now i/ (x,i) and 
we indicate with v (x,f) its actual rate. The symmetric tensor g is the spatial 
metric characterizing the present state of the body; it plays a prominent role 
because in this case the counterpart of 10) of the Lagrangian density is of the 
form 

1 2 

C{x,u,g,v,v,gradv) ^ -p\\v\\ +px{v,v)-pe{g,v,gradv)~pw , (52) 

with some slight abuse of notation. We then find balance equations as conditions 
verifying the relation 

d dr J £ (x, u, g, I/, u,(7radi/)^ =0, (53) 
where S denotes the total variation. 
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To define the variation of the relevant fields, we make use of introduced 
at point 2 of Section 1 and identify 6x with v. We consider a special (though 
wide) subclass of possible vector fields x i — > v (x) characterized by the circum- 
stance that they are purely deformative; in other words, we choose v such that 
skewgradv = 0. 

We then define 

U=o = ivg = 2symgradv = 2gradv, (54) 

dS2 

where f^* means pull back and is thus the autonomous Lie derivative fol- 
lowing the flow V. In analogous way, we put 

5v = Sv + {gradu) v, (55) 



gradSv = gradSu + (gradu) gradv. (56) 
As an intermediate step we notice that 

S e (g, V, gradu) — Se — 
JB Jb 

= J ^2(9ge • gradv + d^e ■ 6u + dgrad^e ■ (^gradSu + (gradu) grad-v^ ^ . (57) 

By developing the variation of (|^ . making use of lt5^ - lt57|) and Gauss 
theorem, we recognize that appropriate balances in the bulk are 



duJO — dxC + div ( 2dgC — (gradu)^ dgmd^jC ) = 0, (58) 



dvC - duC + div (dgraduC) = 0. (59) 
Cauchy stress T is then given by 

T = ~2dgC - (gradu)^ Sa, (60) 

where the actual microstress Sa is defined by 

■Sa = —dgradvC- (61) 

In the case of simple bodies, (|60|) reduces to the well known Doyle-Ericksen 
formula. 

Remark 7. A requirement of invariance of e under the action of SO (3) 
implies that 

skew (2dgC) = e (A^Za + (gradA^) Sa) , (62) 
where Za = —pd^e is the actual self-force and e Ricci's alternating tensor. 
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4.1 Spatial Hamilton equations 



To find appropriate spatial Hamilton equations, we follow the pattern of Section 
1.1. To this end we define spatial canonical standard and substructiiral momenta 
(p and pL respectively) through 

p =du£. , p, ^d-uC. (63) 

Consequently, the spatial Hamiltonian density is given by 

H (x, p, g, iv, fL,gradu) = p • u + /i ■ v—C (x, u, g, u, Vjgradu) (64) 

(with some slight abuse of notation) and has partial derivatives with respect 
to its entries. By evaluating the variation of Ti, taking into account H54|l and 
H56|) . and comparing the result with the variation of £, after making use of the 
balances (|58|l and (|59|l . we obtain the spatial form of the Hamilton equations: 

p = -d^H + div (^2dgH - {gradi>)^ dgradv'Hj , 

u = dpU] (65) 



= -duU + div [dgraduTi) , 
V = dfiU. (66) 

4.2 Spatial Hamilton-Jacobi form 

We may obtain the spatial counterpart of 149(1 by considering a smooth diffeomorphism 

h ■■ (x, p, g, V, Ji,gradi>) i — * (x», p^, g^, i^^, /i^,5radi^,) , (67) 

which generates a new Hamiltonian density 

(x,,p,,g,,iv^,/i^,gradiv,) . (68) 

Now, we may use a generating function S — S {t, x, p*, i/, /i^), and, following 
the same procedure of Section 3, we find that a necessary and sufficient condition 
to assure that 7i, = is 

dtS + n (x, d^S, g, IV, d^S,gradu) = 0. (69) 

4.3 A spatial form of Poisson brackets 

For the spatial Hamiltonian in equations H65|) - (|66|l . taking into account (|54|l - 
1(56(1 . we define a new variational derivative ^ through the relation 

^ (x, p, zv, /i) • V = (^~d^n + div (2dgH - [gradu]^ dgrad^uj^ ■ v, (70) 
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holding p fixed and allowing x to vary, for any v of the kind used in (|54|) - (|56|l . 
Consider a boundary value problem of the type 



(2dgH - [gradu)^ dgradwUj n = d^u (x) , {dgradu'H) n = d^u [v) , on dB, 

(71) 

(where u (x) and u {i>) are the counterparts of the surface potentials U (x) and 
U{v)). 

The total Hamiltonian is now given by H (x, p, /x) = Sb'^ (with some 
slight abuse of notation) and we list only the entries (x, p, Ji) because we 
consider the variational derivative (|65|l below. We consider also arbitrary func- 
tionals F of the type Jg / (x, p, i/, /i), with / a sufficiently smooth scalar density. 

Theorem 3. The canonical Hamilton equation 

F = {F,H}^ (72) 
is equivalent to the Hamiltonian system of balance equations (|65f) - (l66|l with 



OB 



5f SH I SH Sf . 

-7- ■ -r- \dB — 7— ■ -r~ \dB 
ox op ox op 

df sn SH Sf 

B \Su Sfi Sfi Sv 

Sf SH . SH Sf 

\Su Sfi Sf Sfi 

where {•, -j^j is bilinear, skew- symmetric and satisfies Jacobi's identity. 



(73) 



5 Final remarks 

To illustrate possible uses of Theorem 2, we list below some special cases. Anal- 
ogous results accrue from Theorem 3. 

Remark 8. If we choose / = p • v, with v an arbitrary vector, equation 
1)281, ') and the boundary condition l|31() follow immediately from l|36|) . 

Remark 9. Let / = tJ--£,M i'^)' then from ()36|l we get 1)291, ) and the bound- 
ary condition 1)33(1 . 

Remark 10. Let / be of the form 

/ = p-(qx (x-Xo)) + /x-^q, (74) 

with q arbitrary as in previous sections. Consider also, for the sake of simplicity, 
absence of external bulk interactions (the ones accounted for w (x, u)). By using 
(|2H1) and 1231, we obtain from 

e {dpHF'^) ^ A^d^H + {\IA^)' dv^H. (75) 
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These remarks are the HamHtonian counterparts of Remarks 1, 2, 3. Of 
course, Poisson parentheses not only allow one to write in a concise form bal- 
ance equations, but generate articulated geometric structures over the infinite- 
dimensional manifold of mappings showing placements and order parameters, 
and properties of these structures depend also on the geometric properties of 
M. 
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